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摘 E: 为 了 解决 脑 胶 质 瘤 边 界 模糊 、 复 杂 而 导致 的 分 割 不 准确 问题 ,提出 了 一 种 将 灰 度 直方 图 (gray level histogram, 
GLH) 与 改进 的 细胞 自动 机 相 结 合 的 脑 胶 质 瘤 分 割 算 法 。 首 先 ， 对 脑 胶 质 瘤 的 T2 加 权 图 像 和 液体 衰减 反 转 (FLAIR) 


图 像 进行 融合 ; 然后 ， 利 用 友 度 直方 图 特性 增强 脑 胺 质 瘤 区 域 ; 最 后 ， 以 加 权 距 离 为 特征 向 量 用 改进 的 细胞 自动 机 进 

行 分 割 ， 并 得 到 脑 胶 质 瘤 各 组 织 分 割 结 果 。 在 20 组 BraTS2015(brain tumor hes 10 组 临床 脑 胶 
质 瘤 数据 上 进行 分 割 实验 ， 整 个 肿瘤 区 域 及 核心 肿瘤 区 域 的 平均 分 割 准确 率 分 别 达到 90.76% 和 89.73%。 实 验 结 果 表 
明 ， 相 对 于 对 比方 法 ， 所 提 算 法 不 仅 能 更 好 地 分 割 出 对 比 度 明 显 的 胶 质 瘤 区 域 ， 还 一 定 程度 上 解决 了 模糊 胶 质 瘤 区 域 


分 割 不 准确 问题 。 且 该 算法 在 保持 不 增加 算法 复杂 度 的 同时 ， 亦 提高 了 算法 分 割 的 准确 性 和 鲁 棒 性 
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Brain glioma segmentation for multi-modality MR images based on 
gray level histogram and cellular automata 


Yi Fei!, Gong Jing!, Duan Huihong!, Su Guanqun!, Tian Hailong, Nie Shengdong! 
(l. Institute of Medical Imaging Engineering, University of Shanghai for Science & Technology, Shanghai 200093, China; 2. 
Neurosurgery Department, Qilu Hospital of Shandong University, Qingdao Shandong 266035, China) 


Abstract: The fuzzy and complex glioma boundary can cause inaccurate segmentation of the glioma. In order to solve this 
problem, this paper proposed a new glioma segmentation algorithm combining gray level histogram (GLH) with improved 
cellular automaton. Firstly, this method fused T2-weighted and fluid attenuated inversion recovery (FLAIR) MR images of brain 
glioma. Then, it used the histogram feature to enhance glioma region. And, it calculated the weighted distance eigenvector of 
glioma images. Finally, it utilized the improved algorithm of cellular automata to obtain the segmentation result of glioma tissues. 
It separately segmented twenty groups of Brain Tumor Segmentation database data and ten groups of clinical glioma data. The 
average segmentation accuracy rate of the entire tumor area and core tumor area reach to 90.76% and 89.73% respectively. The 
experimental results show that compared with the contrast method, the proposed algorithm can better segment the glioma region 
with obvious contrast. And it solves the problem of inaccurate segmentation due to the fuzzy glioma region to some extent. 
While, it also improves the accuracy and robustness without increasing the complexity. 
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在 外 科 手 术 和 治疗 计划 中 ， 脑 胶 质 瘤 分 割 是 非常 重要 上 且 具有 挑 


战 性 的 一 项 任务 ", 也 是 胶 质 瘤 后 续 评估 的 基础 。 临床 运 / 


中 ， 


脑 胶 质 瘤 是 成 人 中 最 常见 的 原 发 性 脑 肿 瘤 ， 多 由 胶 质 母 医生 通常 根据 临床 经 验 对 每 一 层 肿瘤 图 像 进行 手动 分 割 ， 


这 不 


Bi ETE NG, 且 浸 润 周围 组 织 ”. 脑 胶 质 瘤 占 恶 性 脑 肿瘤 的 81%， ， 仅 耗 费 大 量 的 时 间 精 力 ,而 且 分 割 结果 的 主观 性 强 、 重 复 率 低 。 


包 ; 
发 ; 
年 


授 ， 


元 率 和 死亡 率 较 高 。 胶 质 母 细 胞 瘤 总 体 生 存 率 低 ， 诊 断后 5 ”因而 准确 有 效 地 利用 自动 或 半自动 方法 分 割 脑 胶 质 瘤 就 显得 万 
存活 率 低 至 0.05%~4.7% ”， 平 均 存 活 率 不 超过 14 SA". 为 重要 。 脑 胶 质 瘤 边缘 模糊 内 部 组 织 结构 复杂 , 单一 模 态 的 MR 
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Ds 


像 不 能 提供 足够 的 图 像 信 息 ， 如 图 1 Ata, TIC 图 像 着 重 显 研究 方法 
示 核 心 肿 瘤 区 域 (黄色 箭 头 标志 ) ， 包 含 增 强 和 坏死 部 分 。T2 和 ARARA 
FLAIR 图 像 中 的 高 亮 区 域 补充 显示 水 肿 部 分 (绿色 箭头 标志 )。 1.1 算法 流程 


因此 本 文 利用 TIC、T2、FLAIR 三 种 模 态 的 MR 脑 胶 质 瘤 图 像 本 文 算法 框架 如 图 2 所 示 。 利 用 TIC 图 像 分 割 脑 胶 质 瘤 的 
进行 脑 胶 质 瘤 分 割 。 核心 区 域 , T2 和 FLAIR 的 融合 图 像 分 割 脑 胶 质 瘤 的 水 肿 区 域 。 


本 文 算法 主要 包括 以 下 步骤 : a) 融合 FLAIR 和 T2 两 种 模 态 脑 
Kem A: O 图 像 的 灰 度 直方 图 转换 ，c) 利用 改进 的 细胞 
自动 机 算法 对 脑 胶 质 瘤 进行 分 割 ， 并 得 到 完整 的 分 割 结果 。 


T1C 图 像 


原 多 模 态 MRI 


(a) TIC HA (b) T2 模 态 (c) FLAIR 模 态 
图 1 脑 胶 质 瘤 多 模 态 MR 图 像 


T2 和 FLAIR 图 像 融合 图 像 


细胞 自动 机 (CA) 广泛 应 用 于 图 像 去 噪 ””、 增 强 “””、 分 
割 和 保密 安全 ”等 领域 。 相 比较 其 他 分 割 方法 ， 简 单 易于 实现 
日 可 塑性 高 。 近 几 年 ， 一 些 研 究 者 提出 多 种 利用 CA 算法 进行 


中 瘤 分 割 的 方法 。 如 Desbordes 等 人 "考虑 图 像 的 空间 信息 ， 
引入 26-Moore 型 邻 域 系统 , 并 在 CA 转移 函数 中 添加 邻 域 间 的 
相似 度 信息 ; Sompong 等 人 "提出 FCM-CA 模型 ， 该 模型 利用 改进 细胞 自动 
灰 度 共生 和 矩阵 提取 肿瘤 图 像 特征 ， 然 后 用 FCM 把 这 些 特征 分 
为 四 类 ， 其 中 两 类 特征 代表 肿瘤 区 域 ， 用 这 两 类 特征 对 CA 转 
函数 进行 改进 ; Hamamci 等 人 "提出 CA-Level Set 模型 ， 该 


分 割 结果 i 灰 度 直方 图 转化 图 像 


模型 利用 Grow Cnut 的 转移 函数 并 添加 了 适应 参数 ， 可 以 调整 图 2 算法 流程 图 

转移 函数 到 最 佳 状态 。 但 是 ， 由 于 脑 腕 质 癌 图 像 比较 复杂 上 且 病 1? 图像 融 合 

例 之 间 差 异 较 大 ， 使 参数 优化 困难 ， 进 而 导致 转移 函数 衰减 过 aPC RUE aie Hi le Re ele 
快 的 问题 Sompong 等 人 "提出 邻 域 -加 权 距 离 的 概念 ， 将 图 TAB, SHORTER WEI AB AT (a), IF ER Ue 
像 的 灰 度 信 息 转换 成 加 权 距 离 信息 ， 很 好 地 解决 了 转移 琢 数 变相 生生 和 全 人 全 A 


Aut 


CRER = pa Z 生 过 分 害 | 问 题 。 PEL SIZ. HAS PR 进行 分 割 |， 能 会 造 
减 过 快 的 问题 。 但 是 该 方法 缺乏 邻 域 相似 度 信 息 ， 有 待 改 进 。 生 过 分 割 问题 。 而 单 采用 FLAIR 模 态 图 像 进行 分 割 ， 可 能 会 造 


口 心 


大 部 分 脑 胶 质 瘤 组 织 模糊 不 清 ， 给 分 割 带 来 困难 。 为 了 解决 脑 。 成 分 割 个 完全 的 问题 , 如 图 3 第 一 行 (0) 并 没有 完全 显示 整个 胶 
ORJUSE, Sachdeva 等 人 "提出 结合 图 像 纹理 信 。“ 质 瘤 区 域 ,因此 将 脑 胶 质 瘤 T2 和 FLAIR 两 种 模 态 的 图 像 融 合 ， 
息 的 活动 轮廓 模型 ， 该 模型 利用 肿瘤 的 纹理 信息 和 亮度 信息 来 能 较为 完整 地 显示 出 整个 脑 股 质 瘤 区域 。 并 能 减 噶 单 模 在 股 质 

过 息 的 干扰 。 
变 


— 


E 

建立 静态 场 和 动态 场 ,进而 使 初始 轮廓 准确 地 演化 到 肿瘤 边 ea d 
en ey] LL 2 ee hey eh ys 2H Age 人 (1415 = 1 用 小 波 变化 进行 图 像 融 合 是 现 如 今 的 研究 热点 ””。 小 
但 该 方法 并 未 分 割 出 完整 的 肿瘤 组 织 *Randhewa 等 人 利用 交 。 波 变换 在 提取 图 像 低频 信息 的 同时 ， 还 可 获得 水 平 、 垂 直 和 对 
又 靖 损失 函数 对 肿瘤 边缘 像素 点 进行 准确 的 分 类 。 该 算法 首先 角 三 个 方向 的 高 频 信息 mm， 便 于 对 各 分 解 信息 进行 更 高 层次 的 
用 DNN 对 脑 胶 质 瘤 图 像 进行 初步 的 分 类 ， 然 后 在 交叉 精 损 失 “， 处理， 改善 融合 效果 。 本 文 融 合 的 目的 是 尽 可 能 完全 显示 胶 质 
oa ee CANIN se Sah, 并 保持 清晰 的 边缘 信息 .文献 [17, 18] 实 验 结论 提 到 结 
函数 中 给 予 边缘 像素 点 较 大 的 权重 ， 并 添加 防止 过 分 割 规则 。 。 合 边 缘 算 子 可 以 较 好 地 保留 图 像 的 细节 信息 ， 保 证 图 像 的 清晰 

但 整个 方法 的 计算 量 较 大 ， 在 数据 量 不 足 的 情况 下 ， 会 产生 局 R 
忆 此 本 文采 用 Daubechies 小 波 模型 ， 根 据 Mallat 算法 ,对 


Dk 


最 优 的 问题 。 
两 幅 待 融合 图 像 分 别 进行 分 解 , 得 到 低频 分 量 六 (2 y) 和 高 频 分 
上 述 方法 主要 针对 脑 胶 质 瘤 分 割 中 存在 的 某 个 问题 提出 解 
去 EE mdy), day), diay) 。j 为 分 解 层 数 ，j =3。 h v, 
d 分 


决 方案 ， 且 准确 分 割 出 水 肿 区 域 的 方法 较 少 。 为 了 进一步 提升 AREKE. ae arate WE 
水 肿 区 域 的 分 割 精度 ， 并 解决 转移 函数 衰减 过 快 而 导致 的 欠 分 局 部 方差 融合 规则 对 两 幅 图 像 中 的 相应 分 量 


和 
割 间 题 ， 本 文 提出 了 一 种 结合 次 度 直方 图 特性 和 改进 细胞 自动 图 像 局 部 方差 计算 方法 如 式 (D 所 示 ，m，n 代表 局 部 矩阵 大 
机 的 多 模 态 MRI 脑 胶 质 瘤 分 割 算 法 。 在 BraTS2015 数据 和 临 yy, m=n=3, u ont 值 。 
床 数据 上 的 实验 结果 表明 该 算法 具有 较 好 的 鲁 棒 性 ， 能 够 准确 vor $$ 06,90) 

fey ee FE BOB UI + (1) 
此 方法 在 图 像 局 部 方差 过 小 时 ， 会 使 融合 图 像 失真 。 因 此 


] Canny 算 子 提取 高 频 分 量 的 边缘 ， 再 对 高 频 边缘 分 量 进行 
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局 部 方差 计算 。 
对 于 图 像 A 和 图 像 B 的 低频 分 量 局 部 方差 融合 准则 如 式 
(2) R, FO, y) 是 融合 低频 系数 。 


L(y) Var, (x, y) 2 Var, (x, y) 


FC CA) = fs (x,y) Var, (x, y) < Vary (x, y) 


(2) 


高 频 分 量 融 合 准 则 如 式 (3)-(5) 所 示 , CA y), CEQ) 
分 别 表 示 A、B 图 像 的 高 频 分 量 经 Canny 边缘 检测 后 的 高 频 系 


， 等 : 结合 灰 度 直方 图 和 细胞 自 


示 图 像 位 深 所 能 表示 的 最 大 灰 度 人 
验证 得 出 œ 在 0.006~0.009， 图 


肿瘤 部 分 被 
述 结论 ， 提 取 以 边界 像素 点 ( 红 
Œ) Ch) ras, 可 以 看 出 


IG, j) =cxe 
经 过 以 上 处 理 后 ， 得 到 图 4(g)。 从 图 4(g) Ae 


像 转 换 的 效果 最 好 。 


(—S(i, j)/@) 


Chi naXiv 合 “ERAT 
BAR SIRT 


H, o 是 权重 系数 。 通 过 实验 


凸 出 显示 ， 边 界 更 加 清晰 。 为 了 更 加 直观 的 反应 上 


数 ，VC%(X,y) 和 VCs(X,y) 分 别 是 高 频 边缘 分 量 的 局 部 方差 ; 
FD) (x, y) 为 融合 高 频 系数 ， 包 含 三 个 方向 。 


VC (x,y) = — SS cio D- (3) 


x=] y=1 


VC (x, y) -— FE (Cat D- (4) 


x=l y=1 


(5) 


i a VC4,(x,9) > VC, y) 
FD; (x,y) = j d 4 
Ch (ny) VCE (ny) < VCH (HY) 
最 后 ， 对 得 到 的 融合 系数 进行 小 波 逆 变 换 操 作 ， 得 到 融合 
像 。 如 图 3(c) 所 示 ， 利 用 小 波 融合 模型 可 以 较 适 应 的 将 脑 胶 
质 瘤 T2 和 FLAIR 两 种 模 态 的 图 像 有 效 融合 。 


Ds 


(a) T2 (b)FLAIR (c) 融合 


融合 效果 医 


13 ” 灰 度 直方 图 特性 
TIC 图 像 能 够 较 好 的 显示 脑 胶 质 瘤 中 增强 和 坏死 部 分 ， 但 
级 别 脑 胶 质 瘤 的 增强 或 坏死 区 域 不 明显 ,导致 T1C 图 像 中 的 
区 域 模糊 不 清 。 融 合 图 像 的 水 肿 边 界 和 TI1C 图像 的 模 


肿瘤 区 域 与 正常 组 织 的 灰 度 差 值 较 小 ， 易 导致 分 割 不 准确 。 

度 直方 图 反应 了 图 像 的 灰 度 分 布 情况 与 灰 度 动态 范围 ， 能 够 
用 于 感 兴趣 区 域 的 分 割 与 增强 。 因 此 本 文 利用 脑 胶 质 瘤 图 像 的 

度 直方 图 特性 对 融合 图 像 和 TIC 图 像 的 胶 质 瘤 区 域 进行 增 
强 ， 并 抑制 正常 区 域 。 
图 4 是 以 脑 胶 质 瘤 的 融合 图 像 作 为 输入 图 像 进 行 GLH 
(gray level histogram) 处 理 的 过 程 ， 亦 适用 于 TIC 图 像 。 脑 胶 
质 瘤 区 域 与 正常 组 织 区 域 的 灰 度 存在 差异 ， 且 肿瘤 占 整 个 脑 的 
体积 要 小 于 正常 组 织 。 如 图 4 C) Ce) 所 示 ， 脑 部 正常 区 域 的 
灰 度 值 集中 于 灰 度 直方 图 的 峰 谷 位 置 ， 表 明了 正常 组 织 的 体积 
占 比 较 大 ， 而 肿瘤 和 边界 区 域 的 灰 度 值 集中 于 图 像 的 下 降 及 尾 
部 区 域 。 图 4 Cd) O 直观 地 反映 了 以 上 特性 ， 两 图 中 肿瘤 区 
域 的 灰 度 值 出 现 频率 都 小 于 正常 组 织 的 灰 度 值 出 现 频率 。 利 用 
此 特性 ， 构 造 式 (6) 重新 计算 每 个 像素 点 上 的 值 ，74, J) RER 
像素 点 J) 的 灰 度 值 ，S(i, J) 表示 该 灰 度 值 出 现 的 概率 ， 
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脑 胶 质 瘤 与 了 


色 ) 为 中 心 的 5x5 矩阵 。 如 图 4 
变换 后 图 像 边界 与 背景 区 域 的 灰 度 差 


正常 组 织 的 灰 度 值得 
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(a) 脑 胶 质 瘤 


融合 图 像 的 一 部 分 


(b) 图 (a) 标 记 框 的 灰 度 矩阵 


Hi Danna 
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(e) 原 图 像 灰 度 概 率 图 
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HOMA 


(e) 处 理 后 的 结果 ，om=0.006 
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4 脑 胶 质 瘤 图 像 的 灰 


1.4 改进 的 细胞 自动 机 算法 


度 直方 图 特性 转换 示意 图 


CA 是 基于 像素 的 一 种 算法 


REE? 


演化 分 类 ， 
Hil. 


CA K 


包 ， 在 特定 有 限 的 邻 域 系统 内 根 所 


到 了 有 效 地 抑制 , 相 比 于 原 图 ， 
E 常 组 织 的 边界 更 为 清晰 ， 有 利于 后 续 分 割 。 


， 该 算法 把 图 像 中 的 每 一 个 像 


居 转 移 函 数 进行 


演化 受 周边 细胞 影响 0 站。 且 可 以 实现 图 像 的 多 类 分 


邻 域 像素 点 。 


N(p)={qgeZ":|| p 


—q||.=max | p, =q; Hp 


| Moore 型 八 邻 域 系统 。 如 式 (7) 所 示 , q 是 p 的 


(7) 
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转移 函数 如 式 (8) 所 示 。 


eel, 1,> L 且 1, 是 肿瘤 


8(p,9)=) 1, 
evel, 其 他 


(8) 


上 述 转移 函数 仅 考 虑 两 个 像素 点 之 间 的 灰 度 关系 ， 不 能 准 
确 判 断 边缘 的 位 置 。 如 图 5 所 示 ， 中 心 点 10 与 3x3 邻 域内 像 
素 点 的 灰 度 差 值 相等 ， 无 法 准确 判断 中 心 点 所 属 区 域 。 但 其 邻 
域 像 素 点 所 属 区 域 是 不 同 的 。 图 5 中 ， 以 12 和 8 为 中 心 点 的 
两 个 八 邻 域 是 有 差别 的 ， 且 与 中 心 点 为 10 的 八 邻 域 相似 度 不 
同 。 根 据 以 上 特征 本 文 将 像素 点 之 间 的 关系 转化 为 邻 域 之 间 的 


ji) 


pem 


关系 ， 能 够 有 效 的 区 别 灰 度 相 近 的 像素 点 。 
SSS SR 
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图 5 像素 点 的 空间 信息 

仅 以 灰 度 值 为 特征 向 量 进行 演化 分 类 ， 会 导致 指数 转移 函 
数 衰 减 过 快 。 为 了 降低 转移 函数 的 衰减 ， 并 准确 有 效 地 定位 边 
界 , 本 文 将 采取 将 像素 点 之 间 的 灰 度 关系 转化 为 加 权 距 离 关 系 。 
图 6(b)(c) 分 别 是 以 灰 度 和 均 方差 加 权 距 离 为 特征 向 量 对 (a) 
进行 演化 的 结果 ， 其 中 f=0.05。 可 以 看 出 (b) 中 转移 函数 即使 
在 同一 区 域 也 存在 过 度 衰 减 的 问题 ， 导 致 分 不 清 边界 ， 最 终 衰 
减 结果 为 0.0005; Cc) 解决 了 衰减 过 快 的 缺点 ， 并 能 有 效 地 区 
分 出 边界 和 平坦 区 域 。 
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1 [0.2346 | 0.1003 | 0.045 | 0.2865 


0.8619 | 0.7904 | 0.7223 


0.5488 | 0.6376 | 0.1353 | 0.1003 | 0.3166 


0.9193 | 0.9392 | 0.8251 


0.7760 


0.5382 | 0.7627 | 0.8825 | 0.8167 


0.0008 | 0.0017 | 0.0019 | 0.0019 | 0.0017 | | 0.4682 | 0.4905 | 0.4928 | 0.4979 
[0.4525 |0.4514 | 0.4507 |0.4508 |0.4529 | 
(b) 灰 度 特征 向 量 © 加 权 距 离 特征 向 量 
图 6 转移 函数 衰减 结果 图 
改进 的 细胞 自动 机 算法 是 以 5x5 为 一 个 计算 邻 域 ， 转 化 过 
程 中 同时 考虑 4, 和 9 邻 域 ( 入。) 的 信息 ， 充 分 利用 像素 点 间 
的 空间 信息 。 像 素 点 的 权重 取决 于 该 像素 点 的 邻 域 与 中 心 点 邻 
域 的 相似 程度 ， 本 文 将 邻 域 间 灰 度 差 值 的 离散 程度 作为 像素 点 
之 间 的 相似 度 值 。 演 化 过 程 不 再 仅仅 依赖 于 两 个 像素 点 之 间 的 
灰 度 关系 ， 而 取决 于 两 个 像素 点 之 间 的 加 权 距 离 关 系 。 图 7 是 
计算 像素 点 间 加 权 距 离 的 示意 图 。 
7 (a) 截取 了 实际 图 像 中 的 一 部 分 ， 作 为 转换 的 输入 算 
7 b-c) 计算 入 ,入 ,之 间 的 距离 ad(p,4:) 。 如 式 (9) 
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所 示 ， 其 中 入 , N, 分 别 是 以 p 和 9 为 中 心 点 的 3x3 领域 ， 
GEN,, (x,y) 和 (e, 耻 ) 分 别 表示 p M q; R oreh] 


Ad(p,q, )= IN,-N,| = |[(xtr,,ytr, )-letrn,f+r) (9) 
距离 邻 域 的 均 方 差 能 够 反映 该 邻 域 中 各 距离 值 的 离散 程度 ， 
如 图 7(d) 所 示 , 离散 值 越 小 , 该 邻 域 与 中 心 域 的 相似 度 越 高 。 
R Ad(p.q,) 的 均 方差 ， 如 式 10) 所 示 。 其 中 4 9(p,q) 表示 距 
离 矩 阵 的 平均 值 ，7, 表示 邻 域内 像素 点 数目 ,本 文 取 必 =9 。 


Ad(p,q,) — Ad(p.q,))* 
di (Ad(p 4 (p CONE (10) 


P 


d.(p.4,) = sqrt( 


邻 域 点 9 的 权重 系数 如 式 (11) 所 述 ，Y9 ¢N(p)si=1..n, , 


a 能 够 使 权重 调整 到 一 个 合适 的 位 置 ， 本 文 取 w=0.1 ， 结 果 如 
图 7 Ce) 所 示 。 
et! (Pq) 
WG) == 
OO Festi) (11) 


7 O ~ (h) 表示 利用 式 (9) 和 C11) 将 距离 矩阵 和 权 
重 矩 阵 相 结合 ， 并 求 其 平均 值 即 得 到 相 邻 像素 点 之 间 的 均 方 差 
加 权 距 离 ， 如 式 (12〉 所 示 。 


w, oad (p,q; 
dp (p,q) 2h) (12) 


从 变换 的 结果 可 以 看 出 ， 与 种 子 点 属于 同一 区 域 的 像素 点 
加 权 距 离 较 小 ， 反 之 较 大 ， 这 也 符合 同一 区 域 灰 度 值 相近 ， 不 
同 区 域 灰 度 值 相差 较 大 的 规律 。 但 由 于 像素 点 闻 的 加 权 距离 值 
小 于 灰 度 差 值 ， 因 此 可 有 效 地 降低 转移 函数 的 衰减 程度 ， 所 以 
本 文 将 转移 函数 定义 为 
gpg) =e erel a3) 
均 方 差 加 权 距 离 的 细胞 自动 机 算法 的 流程 如 表 1 所 示 。 
表 1 改进 的 细胞 自动 机 算法 伪 代 码 


Input: GLH transformed image; T1C image; labeled image; 

Method: 

1: for each p el do Il Calculate neighbor weighted distance 
2: for YqeN, do li N pe neighboring cell of p 

3: Calculate the MSE of neighbored distance, d (p,q; ); 

4: endfor 

5: Calculate neighbor weighting of p, W,,(q; ); 

6: Calculate neighbor weighted distance, d pw (PGi 

7 

8 


: endfor 

: for each p € Z do ll Initial seed labeling 
9: if pisaseed, 0% =1; 
10: otherwise, 的 =0; 
11: endif 
12: end for 
13: Do until convergence // Evolution by CA 
14: for each p e7 
15: for Yq, € N, 


16: if g(d (pa; ))x 0g, >0,; 
t+] t, ' 


p t+1 a? t 
18: 9, = gld (Po; ))x6, ; 
19: end if 
20: endfor 
21: endfor 
22: enddo 


Output: segmentation result 
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(Clinical Glioma, CIG) 数据 (7 组 HGG 和 3 组 LGG) 进行 


Input Neighbor feature 

Transformed  p=142 分 割 实 验 。 数 据 库 中 HGG 分 割 结果 包含 水 有 种、 增强 和 坏死 部 
Image q=144;175,1981;102;151,25,28,25 分 , LOG 由 于 增强 和 坏死 部 分 不 易 区 分 ， 因此 分 割 结果 包含 水 
aa [esfas] [esferi 肿 和 核心 (增强 + 坏死 ) 部 分 。 实 际 临床 数据 没有 严格 区 分 增强 
iain. aes 和 坏死 的 金 标准 ， 因 此 亦 只 对 核心 和 水 肿 部 分 进行 分 割 。 
| T” Eee Eek ehe BraTS2015 数据 库 所 有 图 像 都 已 去 掉 颅骨 ， 严 格 与 TI 图 像 对 
一 一 一 a BCE EY [af fiar 齐 并 各 向 同性 的 插值 到 1mm 的 分 辨 率 。 采 集 图 像 的 设备 场 强 
le lal [sels 为 3T， 图 像 大 小 为 240X240 像素 ,并 包含 金 标准 。 同 时 山东 齐 
(a) 原 图 像 输 入 矩阵 (b) pA SR ES 鲁 医院 临床 脑 胶 质 瘤 数据 ， 包 含 资深 影像 科 专 家 手动 分 割 的 金 
Neighbor distance MSE neighbor 标准 , 层 厚 6mm, 设 备 场 强 为 1.5T, 图 像 大 小 为 336X336 像素 。 
distance 实验 中 ， 将 BraTS2015 数据 库 中 的 每 组 数据 随机 选取 10 层 带 
有 脑 胶 质 瘤 的 图 像 进行 分 割 ， 每 组 数据 最 终 的 定量 评价 结果 是 
sp epg pere 10 层 分 割 结果 的 平均 值 ， 山 东 齐 鲁 医院 数据 的 层 厚 比 较 厚 ， 所 

o sansa 


实验 环境 为 mtel i5 处 理 器 ， 主 频 3.2 GHz, 内 存 8 GB. HK 
HEAR (DSC) °°) EDA PER FPR)”, URHE (Sensitivity ) 


> mu 106| 66 | 39 35 | 33 | 30] | 2 | 24 | 54 
99 [117 [123| | s4 [a4] 123] | 81 [114] 126 


22 | 26 | 26 23 | 26 | 26 | | 23 | 26 | 25 (22 作为 定量 评价 指标 。 阳性 (TP), 真 阴性 (TN), 假 阳 性 
人 各 邻 域 与 中 心 域 的 距离 (d) (c) 中 各 矩阵 的 均 方差 CFP) 和 假 阴 性 (FN) 定义 如 下 : TP=RNT 3 TN=RUT ; 
eal FP=RNT ;FN=RNT ,其 中 工 是 胶 质 瘤 实际 位 置 的 像素 总 和 ， 

38 2 24 
13 | 81 | 114| ACs 《 =) = yH. oP gs 
R 是 分 割 结果 的 像素 总 和 。DSC 定义 为 : se ee 假 

20 6 | 16 | 

Neighbor weighting = a 二 阳性 率 定 义 为 : FPR-—__, 敏感 性 定义 为 : Sons. = Tp Fy o 

0.0150 | 0.0101 | 0.0101 [0.0150] 0.0101 | 0.0101, —— 
0.0744 0.6715 | 0.1830 o : 式 (6) 中 的 w 以 0.003 为 步 长 ， 获 取 0.001 至 0.012 的 整 
0.0123] 0.0136 | 0.0101 — 个 肿瘤 分 割 结果 的 DSC 和 FPR 值 ， 如 图 8 (a) 所 示 。 当 w 处 
81 | 114 | 126 | 于 0.006~0.012 If, DSC 值 达 到 最 大 且 趋 于 平稳 , 但 o 在 

23 26 25 


0.009~0.012 时 , FPR 值 上 升 , 这 将 导致 肿瘤 误 分 割 的 概率 增 大 。 


(© 邻 域 权重 © REEE S RRS — HA 因此 @ 取 值 在 0.006~0.009 时 ,可 取得 较 好 的 分 割 结果 。 式 (11) 
= 0.1352 0.1108 10.0000 中 的 a 和 式 (13) 中 的 8 以 0.005 为 起 点 ， 从 0.01 到 0.25 以 


asana faao ta 0.05 为 步 长 ， 分 别 计算 整个 肿瘤 分 割 结果 的 DSC 指标 。 从 图 8 
pagel oma (Cb) 中 可 看 出 ，a 在 0.1 处 取得 最 优 值 ， 而 8 对 结果 的 影响 较 
03004 | 0.0604 | 0.16 Istance sega < = ab plete popo 
PAT 小 ,这 也 证 明了 采用 加 权 距 离 作 为 特征 向 量具 有 较 好 的 鲁 棒 性 。 
1 


TARL be a 1.2463 | 3.8440 | 4.0871 09 
0.8 
: 3.3035| 0 | 1.5284 07 
0.5 
— 12.4415] 11.9295] 11.9263 0s ae 


0.3 = FPR 
6.0266 | 76.5486 23.0579 0.2 
0.2829 | 0.3534 | 0.2517 号 
0.001 0.003 0.006 0.009 0.012 
(g) DA ER (h) 邻 域 加 权 距 离 (a) 参数 © 的 DSC 和 FPR 指标 
图 7 均 方差 加 权 距 离 计 算 流 程 0.95 5 
0.9 i 
2 ”实验 结果 与 分 析 rey 
R tE zZ 0.75 4 
2.1 算法 的 准确 性 分 析 | = 
为 了 验证 所 提出 算法 的 有 效 性 , 本 文采 用 20 组 BraTS2015 0.65 ] =i 
数据 库 中 的 临床 数据 (包含 10 组 高 级 别 脑 胶 质 瘤 (High Grade or 0.01 0.05 0.1 015 02 0.25 
Glioma, HGG) 数 据 和 10 组 低级 别 脑 胶 质 瘤 (Low Grade Glioma, (b) 参数 a.B 的 DSC 指标 
LGG) 数据 ) 以 及 山东 齐鲁 医院 提供 的 10 组 实际 临床 脑 胶 质 瘤 8 参数 优化 
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(a) GLH 融合 


图 9 HGG_0009、HGG_147 分 割 结果 示意 


(d) 核 心 区 域 


R] 


(a) GLH we 
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10 LGG_0015, LGG_470 分 割 结果 示 意 


(d) 核 心 区 域 


Ka 


(a) HGG_0009 (b) HGG 147 (c) LGG 0015 (d) LGG 470 


PS 


像 上 分 割 核心 肿瘤 部 分 。 高 级 别 脑 胶 质 瘤 核心 
水 肿 区 域 模糊 不 清 。 对 融合 图 


9 显示 的 是 两 侦 


在 GLH 转换 后 的 融合 


分 割 结果 与 金 标准 对 比 民 
高 级 别 脑 胶 质 瘤 病 例 和 | 结果， 分 别 
到 像 上 分 割 整体 肿瘤 部 分 ， 在 TIC 原 图 


区 域 比较 明显 ， 


质 瘤 区 域 , 有 效 抑 


E 常 组 织 。 


像 进行 GLH 254 
图 10 是 两 例 低级 


AY E A h h 


例 的 分 者 


1 结果。 对 于 低级 别 脑 


体 区 域 和 核心 肿瘤 
RUT Aaa 


区 域 影像 特征 都 不 明显 ， 无 法 状 


别 脑 胶 质 瘤 病 
区 中 
KI 


因此 对 融合 


像 和 TI1C 图 像 都 进行 灰 度 直方 图 特 1 


0 


然后 分 别 在 图 10(a)(b) 上 分 


又 域 ， 如 图 10(c)(d) 所 示 。 


1 出 低级 别 朋 
图 11 第 一 行 是 将 各 ; 


或 还 是 水 


(a)(b) 所 示 。 


疝 胺 质 瘤 的 整体 和 核心 
府 例 脑 胶 质 瘤 整 


强 区 域 ， 蓝 色 代表 # 


的 分 割 结果 按照 
。 前 两 列 粉色 代表 水 肿 区 
区 域 。 后 两 列 粉色 代表 水 肿 区 域 ， 绿 色 


崭 胶 质 瘤 的 空间 位 
或 ， 绿 色 代表 增 


ji 
ok 


区 域 。 第 二 行 是 各 病例 的 金 标准 ， 可 看 出 各 例 的 分 割 


标准 


为 了 进 
HE AE h 
55 CA-L 


人 


Cut 是 CA 


- SJA HOT) 
Chi c HHT] 


相似 度 较 高 ， 达 到 了 较 准 确 分 割 的 目的 。 
步 证 明 本 文 算法 的 有 效 性 ， 本 文 分 别 对 30 组 数 
瘤 分 割 结果 与 核心 肿瘤 分 割 结果 进行 了 定量 评价 ， 


evel set 模型 中 和 Grow Cut 模型 进行 对 比分 析 。Grow 


的 经 典 算法 , CA-Level set 模型 将 改进 的 细胞 自动 机 


模型 相 结 合 ， 组 成 了 一 种 有 效 的 分 割 方法 。 

从 图 12 可 看 出 , 本 文 算法 取得 了 较 好 的 分 割 结果 。 数据 库 
中 的 高 级 别 和 低级 别 病例 的 整个 肿瘤 分 割 结果 的 DSC 平均 值 
分 别 为 0.93 和 0.92。 高 级 别 脑 胶 质 瘤 的 核心 肿瘤 部 分 和 增强 部 


割 结果 


的 DSC 平均 值 分 别 达 到 了 0.95 和 0.90。 而 低级 别 


心 肿瘤 分 割 结果 的 DSC 平均 值 为 0.89。 由 于 高 级 别 脑 胶 质 
瘤 的 增强 和 坏死 部 分 较 明显 ， 而 低级 别 脑 胶 质 瘤 特征 不 明显 ， 
各 部 分 组 织 都 比较 模糊 。 所 以 高 级 别 脑 胶 质 瘤 的 分 割 结果 相对 


分 胶 质 瘤 4 


DSC 


1 
3 
a 09 


0.8 
0.7 


0.6 
0.5 


0.4 
0.3 
0.2 


0.1 
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较 好 。 实 际 临 床 中 的 脑 胶 质 瘤 没 有 去 除 脑壳 且 病 例 较 复 杂 ， 部 
长 于 脑壳 边缘 ， 因 此 较 数 据 库 中 病例 的 分 割 结果 
H, 整个 脑 胶 质 瘤 区 域 和 核心 区 域 分 割 结果 的 DSC 平均 
值 分 别 为 0.88 和 0.86. 
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(a) BraTS2015_ 高 级 别 脑 胶 质 瘤 分 割 结果 的 DSC 指标 
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(b) BraTS2015_ 低 级 别 脑 胶 质 瘤 分 割 结果 的 DSC 指标 
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(c) 实际 临床 脑 胶 质 瘤 分 割 结果 的 DSC 指标 
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表 2 将 本 文 算法 与 CA-Level set 模型 和 Grow Cut 模型 进 于 对 比方 法 。 其 中 低级 别 脑 胶 质 瘤 分 割 结果 的 平均 DSC Ail Sens. 
行 了 定量 对 比 。 从 表 2 可 发 现 本 文 算法 要 明显 优 于 其 他 两 种 贷 值 分 别 达到 了 0.902 和 0.89， 实 际 临 床 脑 胶 质 瘤 分 割 结果 的 平 
法 .对 于 数据 库 中 高 级 别 的 脑 胶 质 瘤 , 因 其 组 织 差 异 比较 明显 ， 均 DSC 和 Sens. 值 分 别 为 0.869 和 0.852。 同 样 ， 对 比方 法 对 
所 以 整体 部 分 和 核心 部 分 的 分 割 结果 都 比较 好 。 本 文 算法 的 LGG 和 CIG 的 整体 肿瘤 部 分 和 核心 肿瘤 部 分 的 分 割 结果 值 比 
DSC 和 Sens. 的 平均 值 分 别 达 到 0.937 和 0.916， 并 且 各 病例 之 HGG 略 低 。 且 对 CIG 的 分 割 结果 值 最 低 。 

间 的 差距 和 假 阳性 率 都 要 小 于 其 他 两 种 方法 。 由 于 核心 肿瘤 藤 CA-Level set 模型 可 以 调节 转移 函数 达到 一 个 较 好 的 分 割 
边界 较为 清晰 , 因此 对 比方 法 对 整个 肿瘤 分 割 结果 的 DSC 值 要 状态， 而且 利 用 水 平 集 算 法 对 结果 进行 优化 ， 因 此 分 割 结 果 比 
小 于 核心 肿瘤 的 DSC 值 . 对 于 数据 库 中 低级 别 脑 胶 质 瘤 和 临床 Grow Cut 模型 的 分 割 结果 好 。 此 外 ，Grow Cut 模型 存在 过 早 售 
脑 胶 质 瘤 ， 由 于 低级 别 脑 胶 质 瘤 的 水 肿 部 分 与 核心 肿瘤 部 分 的 ” 止 演化 的 问题 ， 会 导致 从 分 割 。 且 CA-Level set 模型 采用 水 平 
边界 较为 模糊 ， 且 实际 临床 的 脑 胶 质 瘤 比 较 复杂 ， 图 像 质量 不 集 平 滑 边缘 ， 该 方法 对 初始 轮廓 要 求 较 高 ， 且 参数 优化 较 难 ， 
高 ， 无 法 正确 判断 边界 。 因 此 本 文 算法 对 低级 别 脑 胶 质 瘤 和 临 。” 易 导致 边缘 定位 不 准确 。 上 述 问 题 会 直接 导致 DSC 和 Sens. 值 
床 脑 胶 质 瘤 分 割 结果 的 平均 DSC 和 Sens. 值 略微 低 于 数据 库 中 ”降低 。 


jy, 
的 高 级 别 胶 质 瘤 ，FPR 值 也 略微 提高 。 但 分 割 结果 都 要 明显 优 
表 2 三 种 算法 对 HG、LG 和 Cl 的 整体 和 核心 肿瘤 部 分 分 割 结果 的 定量 评价 指标 (平均 值 土 均 方 差 ) 


方法 Whole Core 
DSC FPR Sens. DSC FPR Sens. 

本 文 算法 0.927+0.03 0.028+0.01 0.909+0.0 0.94740.0 0.032+0.02 0.923+0.04 
HGG CA-Level set 0.705+0.07 0.090+0.09 0.690+0.1 0.786+0.0 0.091+0.08 0.749+0.05 

Grow Cut 0.687+0.10 0.164+0.03 0.682+0.1 0.762+0.1 0.120+0.09 0.696+0.11 

本 文 算法 0.916+0.03 0.039+0.02 0.905+0.0 0.88740.0 0.076+0.05 0.875+0.06 
LGG CA-Level set 0.72140.11 0.046+0.07 0.707+0.1 0.705+0.1 0.04740.07 0.685+0.16 

Grow Cut 0.633+0.13 0.010+0.02 0.629+0.1 0.655+0.1 0.058+0.06 0.653+0.16 

本 文 算法 0.880+0.01 0.040+0.03 0.851+0.0 0.858+0.0 0.029+0.02 0.852+0.03 
CIG CA-Level set 0.645+0.13 0.120+0.03 0.646+0.1 0.610+0.0 0.039+0.07 0.633+0.08 

Grow Cut 0.578+0.15 0.016+0.05 0.594+0.1 0.607+0.0 0.079+0.14 0.622+0.06 


如 图 13 所 示 ， 本 文 算 法 的 DSC 值 均 高 于 其 他 两 种 算法 ， 且 在 某 些 病例 中 的 FPR 值 低 于 本 文 算法 。 因此 通过 分 析 可 得 出 ， 
FPR 值 较 低 且 更 稳定 。 由 于 高 级 别 脑 胶 质 瘤 结 构 比 较 复 杂 ， 通 CA-Level set 模型 的 分 割 结果 要 优 于 Grow Cut 模型 。 而 本 文 算 
常会 与 脑室 相连 接 ， 导 致 对 比方 法 分 割 结果 的 FPR 值 较 高 。 对 法 可 进一步 解决 CA-Level set 模型 与 Grow Cut 模型 分 割 过 程 
于 低级 别 和 实际 临床 的 脑 胶 质 瘤 ， 胶 质 瘤 的 对 比 度 较 低 且 临床 ” 中 存在 的 问题 ， 并 取得 较 好 的 分 割 结果 。 从 以 上 实验 结果 可 看 
胶 质 瘤 的 病症 表现 多 样 且 复杂 ， 对 比方 法 易 产 生 从 分 割 结果 ， 出 ， 本 文 算法 能 够 较 准确 地 分 割 脑 胶 质 瘤 ， 并 且 和 鲁 棒 性 较 高 。 
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一 一 本 文 算法 “一 "一 CA- Levelset —*— Grow Cut 


图 13 整体 脑 胶 质 瘤 分 割 结 果 的 DSC 和 PPR 指标 对 比 图 
本 文 算法 的 鲁 棒 性 主要 体现 在 对 于 参数 变化 的 鲁 棒 性 和 对 ” 约 在 0.5 左右 会 得 到 最 好 的 分 割 结 果 。 

于 病例 多 样 化 的 鲁 棒 性 。B 是 三 种 算法 的 核心 参数 ， 用 来 控制 统计 各 算 法 对 三 种 病例 分 割 结果 重合 率 的 均 方差 ， 均 方差 

转移 函数 演化 趋势 ， 影 响 分 割 结果 。 因 此 记录 参数 B 在 0.1 至 可 以 反映 出 每 个 分 割 结果 之 间 的 差距 ， 数 值 越 大 ， 差 距 越 大 ， 

0.7 之 间 变 化 时 ， 各 算法 分 割 结果 的 DSC 值 ， 如 图 14(a)。 可 看 ”和 鲁 棒 性 越 差 。 图 14(b) 中 ,从 纵向 看 本 文 算法 分 割 结果 的 均 方 差 

出 本 文 算 法 在 保持 较 高 重合 率 的 同时 ， 分 割 结果 也 基本 不 受 参 。 要 小 于 其 他 两 种 算法 ，Grow Cut 分 割 结果 的 差距 最 大 ， 在 单 类 

数 B 的 影响 ， 线 形 较 平稳 。 而 对 比方 法 随 p 值 变化 而 变化 ， 大 ”病例 中 本 文 算法 的 鲁 棒 性 好 于 对 比方 法 。 从 横向 看 对 于 不 同类 
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的 病例 ， 本 文 算法 分 割 结果 的 均 方差 相差 较 小 ， 


表明 本 文 算法 


在 对 于 病例 多 样 性 方面 鲁 棒 性 较 好 。 而 其 他 两 种 算法 对 各 类 病 


下 一 本 文 算法 一 -CA-Levelset ——Grow Cut 


(a) 参数 B 


2.2 算法 复杂 度 分 析 
2.2.1 运算 复杂 度 分 析 


对 各 算法 结果 的 影响 


本 文 算法 主要 


三 个 部 分 组 成 ， 因 而 运算 复杂 度 主要 集中 


即 OU) ;图 像 灰 度 直 方 图 转化 部 分 复杂 


在 这 三 个 步骤 上 。 图 像 融合 过 程 中 在 矩阵 融合 部 分 会 产生 一 个 
O(n?) 的 复杂 度 , 如 式 (2) 和 (5) 所 示 , 其 他 行 复 杂 度 为 常数 ， 


FEN OW’); 改进 的 CA 


算法 中 计算 均 方差 加 权 距 离 和 初始 化 种 子 模板 这 两 步 的 运算 复 
ALAA OW) ,而 CA 演化 部 分 的 运算 复杂 度 为 002) 。 因 此 本 


文 算法 的 运算 复杂 度 为 002) o 
该 方法 没有 图 像 预 处 
分 的 运算 复杂 度 为 ON) 。 
法 的 复杂 度 为 0(n*) ， 水 平 集 演化 部 分 
此 整个 模型 的 复杂 程度 为 O(n”) o 


Grow Cut 源 于 CA 算法 ， 虽 然 
的 步骤 ， 但 是 在 最 后 像素 点 演化 分 类 部 
CA-Level set 模型 的 前 半 部 分 CA 算 


的 复杂 度 亦 为 002) ， 因 


三 种 算法 的 运算 复杂 度 没 有 


差别 ,但 是 运算 时 间 有 所 差别 , 如 表 3 第 三 列 所 示 。 经 典 的 Grow 


Cut 算法 因 其 对 整 幅 图 像 进行 处 理 ， 
趣 区 域 ， 导 致 运行 时 间 过 长 。 


耗费 大 量 的 时 间 在 非 感 兴 
本 文 算 法 和 CA-Level set 算法 只 


对 感 兴趣 区 域 进行 处 理 ， 因 此 算法 
文 算法 运 


Ol 


an CA § 法 会 耗费 


的 运行 时 间 减 少 了 很 多 。 本 
行 时 间 主 要 耗费 在 改进 的 CA 分 割 步 又 上 ， 因 为 改进 
后 的 CA 需要 将 图 像 的 灰 度 信息 转化 为 加 权 距 离 信 息 ， 相 对 于 
定 的 时 间 。 而 CA-Level set 算法 中 水 平 集 


衣 8, 等 


合 灰 度 直 方 图 和 细胞 自 蕊 


例 的 分 割 


吉 果 差别 较 大 ， 尤 其 对 于 低级 别 脑 胶 质 


瘤 ， 说 明 对 比 


方法 对 低级 别 脑 胶 质 瘤 分 割 的 鲁 棒 性 较 差 。 
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O) 不 同类 胶 质 瘤 分 割 结 果 的 差异 比较 


图 14 各 算法 鲁 棒 性 比较 
演化 至 收敛 时 需要 耗费 时 间 ， 


平均 运算 时 间 比 本 文 算法 稍 短 。 


实验 中 发 现 ， 本 文 算 法 和 CA-Level set 算法 的 运行 时 间 与 


ROI 区 域 大 小 有 关 ， 
个 值 作 比 较 ， 如 图 15 (a) 所 示 。 


寻 此 取 实 验 中 最 大 和 最 小 ROT 之 间 的 五 
可 以 看 出 两 种 算法 随 ROI 的 


增 大 运行 时 间 也 增 大 。 由 于 本 文 算法 的 运算 量 主要 集中 在 像素 
点 的 分 类 上 ， 因 此 本 文 算法 的 运行 时 间 易 受 ROI 大 小 的 影 
响 ， 图 15 Cad 可 以 看 出 当 输 入 ROT 很 小 的 时 候 运行 时 间 较 


短 ， 输 入 越 大 运行 时 间 增 长 越 快 。 
Level set 部 分 与 ROI 大 小 联系 不 大 ， 


而 CA-Level set 算法 中 
主要 联系 在 CA 部 分 ， 随 


着 ROI 的 增 大 而 增 大 ， 比 本 文 算 法 增 大 的 速度 慢 。 


2.2.2 空间 复杂 度 


度 分 析 


三 种 方法 的 空间 复杂 度 相 当 ， 如 表 3 最 后 一 列 。 本 文 算 法 
在 运行 的 过 程 中 能 够 及 时 的 释放 空 


存储 量 。 


间 ， 不 储存 过 程 结果 ， 减 少 


相对 于 CA-Level set 算法 平均 内 存 需求 稍 低 。Grow 
Cnut 算法 相对 简单 ， 
数据 相对 较 低 ， 对 内 存 需 求 最 低 。 同 样 ， 
区 域 大 小 做 了 比较 ， 如 图 15 (b) 


空间 复杂 度 不 高 ， 运 行 过 程 中 产生 的 存储 
对 内 存 消耗 与 ROI 
所 示 。 随 着 ROI 的 增 大 内 存 


需求 会 增加 ， 但 是 到 了 一 定 的 数值 ， 该 变化 趋 于 平缓 。 总 的 来 


说 ， 三 种 算法 对 内 存 的 需求 都 不 高 。 


表 3 三 种 算法 复杂 度 的 对 比 


运算 复杂 度 平均 运算 时 间 /s 平均 内 存 用 量 /MB 
本 文 算法 O(n’) 98.85 23.84 
CA-Level set O(n’) 93.32 24.20 
Grow Cut O(n’) 212.75 20.96 
120 28 4 
110 26 4 
100 + 24 
90 Sz 
g 2 204 
= so Z 18 
B ng 让 
g 164 
60 ] 一 本文 算 法 14 4 到 -本 文 算法 
50 4 —e—CA-Level set 124 —e—CA-Level set 
40 + + + T 5 10 
43.61 54,52 65.42 72.69 83.59 43.61 54.52 65.42 72.69 83.59 
ROI Side/pixel ROI Side/pixel 
(a) 时 间 变化 (b) 内 存 变化 
图 15 ”时间 和 空间 复杂 度 随 ROT 区 域 变化 区 
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本 文 提 出 的 灰 度 直方 图 结合 细胞 自动 机 的 多 模 态 MRI 脑 
BCA SEI, PG, RF T2 Al FLAIR 两 种 模 态 的 脑 胶 质 瘤 
像 相 融合 ， 可 以 补充 显示 脑 胶 质 瘤 水 肿 部 分 且 可 以 抑制 周围 
和 元 余 信 息 。 然 后 ， 采 用 脑 胶 质 瘤 图 像 的 灰 度 直方 图 特性 
瘤 部 分 凸 出 显示 ， 一 定 程度 上 解决 了 因 肿 瘤 区 域 模糊 而 导致 分 
割 不 准确 的 问题 。 其 次 ， 将 图 像 的 灰 度 信息 转换 为 加 权 距 离 
息 ， 减 绥 了 转移 函数 的 衰减 速度 ， 从 而 使 本 算法 能 有 效 
边界 和 平坦 区 域 。 并 且 充 分 利用 了 像素 点 的 空间 信息 ， 能 够 更 
好 地 分 辨 相似 区 域 。 通 过 实验 结果 证 明 ， 本 文 算法 能 够 较 准确 
地 分 割 脑 胶 质 瘤 ， 并 且 更 加 稳定 。 
站 胶 质 瘤 分 割 是 胶 质 瘤 放射 治疗 、 临 床 手 术 规 划 和 良 恶性 
评估 的 基础 。 在 临床 上 ， 医 生 主 要 依靠 手术 导航 、 荧 光 剂 染色 
等 辅助 技术 来 确定 脑 胶 质 瘤 的 大 概 位 置 。 但 是 ， 现 存 技术 不 成 
熟 ， 对 于 一 些 低级 别 或 复杂 的 脑 胶 质 瘤 ， 胶 质 瘤 区 域 与 正常 组 
织 影 像 特 征 相差 不 明显 ， 不 能 准确 勾画 胶 质 瘤 轮 廊 ， 无 法 满足 
需求 。 这 也 是 现 阶段 研究 者 急需 解决 的 一 个 问题 。 本 文采 
像 的 灰 度 直方 图 特性 凸显 模糊 肿瘤 区 域 ， 并 抑 
组 织 ， 一 定 程度 上 解决 了 该 问题 。 虽 然 两 种 模 态 的 图 像 融 合 可 
以 减少 胶 质 瘤 周 围 的 元 余 信 息 ， 但 是 ， 灰 度 直 方 图 特性 转换 后 
的 图 像 中 脑 胶 质 瘤 边缘 部 分 还 会 存在 小 部 分 干扰 信息 ， 影 响 分 
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割 准确 率 。 如 何 减弱 这 些 干扰 信息 ， 并 最 大 程度 的 凸显 模糊 胶 
质 瘤 区 域 是 接 下 来 研究 的 重点 。 算 法 的 自动 化 程度 在 临床 上 也 
很 重要 ， 自 动 化 程度 越 高 ， 就 可 以 更 多 的 减少 人 为 干预 ， 避 免 


分 割 结果 受 医 生 主 观 性 影响 ， 减 少 误 
ARER. 


诊 率 。 因 此 后 续 研 究 应 结 


合 临 床 应 用 ， 
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